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Abstract 

We present a theory of quantum circuits based on logical qubits encoded in chirality of electron 
spin complexes in lateral gated semiconductor triple quantum dot molecules with one electron spin 
in each dot. Using microscopic Hamiltonian we show how to initialize, coherently control and 
measure the quantum state of a chirality based coded qubit using static in-plane magnetic field 
and voltage tuning of individual dots. The microscopic model of two interacting coded qubits is 
established and mapped to an Ising Hamiltonian, resulting in conditional two-qubit phase gate. 
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I. INTRODUCTION 



There is currently interest in exploiting electron spin for nano-spintronioi and quantum 
information processing^-. This is partly motivated by electron spin long coherence times^ 
and availability of scalable semiconductor technology. In the simplest approach, a physical 
qubit is identified with the two states of an electron spin, which can be manipulated by 
applying local magnetic fields. Much progress has been achieved using micro-magnet^ 
technology with electron spin qubits. An alternative approach is to encode a logical qubit 
in a two level system constructed with spin complexes 3 . This includes singlet-triplet two 
electron qubit£ and a logical qubit encoded in a degenerate ground state of a three spin 
complex—. As suggested by DiVincenzo et al. 10 , in the framework of a Heisenberg model 
for a logical qubit encoded in a linear chain of three spins, such a logical qubit can be 
manipulated by the control of the exchange interaction between pairs of spins. The ability 
to manipulate spin state with a voltage is related to the relation between the orbital and spin 
part of the many-electron wavefunction. A detailed microscopic model of a three electron 
complex in a single semiconductor quantum dot (QD), exploring the orbital and spin relation, 
has been investigated and compared with experiment by one of us 11 . A similar microscopic 
model of three electrons localized in a lateral gated triangular triple quantum dot molecule 
in the plane of a GaAs/GaAlAs heterojunction was proposed and explored by one of us in 



Ref. 
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13 . The coded qubit was identified with chirality of three electron spin complex, or 



equivalently, two possible directions of a minority spin motion. Other proposals to implement 
coded qubits with both triple quantum dots^ (TQDs) and atom traps^^also exist. 

The advantages of working with a TQD-based coded qubit are two- fold. First, every 
quantum gate can be implemented electrically. In such a scheme, magnetic field will only 
be used for initialization and (or) measurement of a coded qubit. Second, the coded qubit 
involves Decoherence Free Subspac o M ^ £~— (DFS) and is immune to channels of collective de- 
coherence. This reduces decoherence of a TQD-based coded qubit due to charge fluctuations 



as discussed in Ref. 
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Other interesting phenomena involving TQD molecules include non- Fermi liquid behavior 
when coupled to leads^i^, potential for generating maximally entangled three-partite GHZ 
and W states^^and manipulation of total spin^i. 

In this work we develop a microscopic theory of quantum circuits based on coded qubits 
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encoded in chirality of electron spin complexes in TQD. We use a combination of Linear 
Combination of Harmonic Orbitals - Configuration Interactions (LCHO-CI) 13 , Hubbard and 
Heisenberg models to determine a set of optimal conditions for single qubit operations and 
describe the two qubit gate. We show that there exists an in-plane magnetic field direction 
and magnitude optimal for single qubit operations. However, the magnetic field will rotate 
qubit states in an undesirable way; and we show how this rotation can be controlled by 
tuning voltages on the gates. Once the optimal sets of magnetic fields and voltages have 
been obtained, the exact diagonalization techniques are used to establish and verify an 
effective two-qubit Hamiltonian. It is shown that the effective two qubit interaction is Ising- 
like, leading to a two-qubit phase gate. The present work establishes both single qubit 
and two-qubit operations necessary for performing quantum computation using TQD-based 
coded qubits. 

The paper is organized as follows. Sec. fl] contains introduction. Sec. [TT1 describes 
the quantum circuit based on chirality of 3 electron spin complexes and computational 
methodology. In Sec. IHIl we present a definition of a TQD-based coded qubit, discuss 
initialization, single qubit operations using voltages, and measurement of the state of the 
coded qubit. In Sec. HVt we present an effective Hamiltonian of 2 coupled coded qubits 
and show that it can be canonically transformed to an Ising interaction. Sec. |V] contains 
summary. 

II. THE MODEL 

Fig. (pji) shows schematically layout of quantum circuits based on coded qubits encoded in 
chirality of electron spin complexes in triangular TQDs. The circles denote individual lateral 
quantum dots formed in the 2D electron gas (2DEG) at the hetero junctions of AlGaAs/GaAs 
by metallic gates on the AlGaAs surface. The gates are set to confine a single electron in each 
dot denoted by an arrow. Additional gates (not shown here) are used to control tunneling 
between dots in the same TQD molecule, shown schematically as solid lines. Dashed lines 
indicate tunneling between neighboring TQDs. The tunneling between any pair of QDs is 
responsible for exchange interaction of electron spins localized on each QD. The brackets 
indicate two TQDs isolated from the rest of the circuit. It is assumed that any number 
of TQDs can be isolated from the rest of the circuit by turning off the boundary exchange 
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interactions with gate voltages. 

For comparison, Fig.(J]l)) shows the circuit composed of linear, instead of triangular, 
TQDs, similar to the linear chain of spins proposed by DiVincenzo et al., in which the coded 
qubits in the chain are implemented with only two exchange interactions. 

Fig. (JTb) is another possible architecture studied by Weinstein^ et. al. for spin system. 
However, in this design, bringing TQDs close together induces interactions between quantum 
dots beyond the ones indicated by dashed line. 

In this work we focus mainly on triangular TQDs-based quantum circuit shown in 
Fig. ([Ik). Since electrons are well localized in each QD, a system of two TQDs in a chain 
can be very well described by an extended Hubbard model™. With cj a (ci a ) electron creation 
(annihilation) operator for the electron with spin a = ±1 on the i — th QD, the Hubbard 
Hamiltonian reads: 

2 3n 2 3n 

n=l i=3n— 2 a n=l i,j=3n—2 a 

. 2 3n 2 3n 

n=l i,j=3n—2 n=l i=3n— 2 

+ (4c 4 a + 4/3.) + V'p 3 p 4 , (1) 

a 

where n = 1,2 labels the two TQD molecules from left to right, indices i,j range from 
3n — 2 to 3n label each QD from left to right in the n — th TQD molecule. The intra- 
TQD Hubbard parameters tij, V, U and E ia are the tunneling matrix element between 
the i — th and the j — th QDs, the off-site Coulomb interaction parameter between any 
two QDs in the same TQD molecule, the on-site Coulomb interaction strength for any 
QD, and the on-site energy for the i — th QD respectively. The on-site energy depends on 
spin and magnetic field which is applied in the plane of a TQD molecule. The Hubbard 
parameters t' and V represent the inter-molecular tunneling matrix element and inter- 
molecular Coulomb interactions. h ia = c\ a Ci a is the number operator for the a electron on 
the i — th QD. pi = h,^ + is the electron charge density operator on the i — th QD. 
When we discuss triangular resonant TQDs, we drop the subscripts on all the parameters. 
In our model, we consider the parameters corresponding to regime of strong correlations: 
4 < t u « V < V < U. 
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The intra-TQD Hubbard parameters, Ei, Uj, V, and U, in Eq.flTJ are obtained from a 
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microscopic calculation for single TQD based on LCHO-CI method as explained in Ref. 
The inter-TQD Hubbard parameter V is taken to be the direct Coulomb interaction between 
2 charges localized on adjacent edge dots of two neighboring TQD molecules. 

The energy spectrum and eigenstates of the Hubbard Hamiltonian for one and two 
TQD molecules are obtained using Configuration Interaction technique. For a given num- 
ber of electons N e we construct all possible configurations \k = iN e CN e , ■-, £202, iitfi) = 
c\ N aN ■■■c\ 2a2 c\ iai |0), build Hamiltonian matrix in the space of configurations, and diago- 
nalize it numerically. 

At half-filling, the low-energy spectrum of the Hubbard model can be approximated by 
a spectrum of a Heisenberg model2&2&22 describing electron spins localized in each dot: 

Hheis = ^ JijSi ' Sj. (2) 

i<j 

The exchange interactions, J^-, for the TQDs can be expressed in terms of tunneling matrix 
elements and quantum dot energies: 

Jij = 2lUjl2 {u-V + iEi-Ej) + U-V-iEi-Ej)) ■ (3) 
The exchange interaction can be controlled by either tuning the tunneling matrix element 

tij by, for example, additional gates controlling the height of the tunneling barrier, or by 

biasing the dots and changing their on-site energy E^y 



III. CODED QUBIT ENCODED IN CHIRALITY OF ELECTRON SPIN COM- 
PLEX IN A TRIPLE QUANTUM DOT MOLECULE 

In this section, we discuss a single coded qubit shown in Fig.([2h): its preparation, ini- 
tialization, operation and measurement in the presence of a lateral magnetic field in the y 
direction. 

The qubit is encoded in quantum states of a three electron spin complex in a fully 
symmetric and half-filled TQD. An example of one of the three possible configurations 
|t44) = cLcLcgjO > with Sy = —1/2 is shown schematically in Fig.(j2k). The three con- 
figurations with the minority spin on QD 1, 2 or 3 form a doubly degenerate ground state 
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with total spin S = 1/2 and an excited state with total spin S = 3/2 separated from the 
ground state by J e ff, as shown in Fig.(j2b). 

The doubly degenerate ground state with fixed S y = —1/2 forms an effective two level 
system. We identify the coded qubit levels \q + ) and \q_) with 



where s = ± and \n) = The three spin complex is characterized by chirality 

X = (Si x S2) - S3 which measures the degree of collinearity of the three spins. The two coded 
qubit levels, Eq.(@|, are the eigenstates of chirality operator with eigenvalues x — ±v3/4. 
These two states also can be more intuitively characterized by minority spin moving either 
to the left or to the right as in Resonant Valence Bond (RVB) plaquette^"— . 

We note that in the absence of the magnetic field the coded qubit states in S y — 1/2 
subspace are degenerate with coded qubit states in S y = —1/2 subspace. If we are to work 
in the computational space corresponding to S y = —1/2, any process which flips the electron 
spin will remove the coded qubit from its computational space. In order to separate the 
computational Hilbert space S y = —1/2 from S y = +1/2 subspace, a magnetic field B y is 
applied along the y direction as shown in Fig(T2h). We avoid applying magnetic field applied 
in the z direction, because it activates undesirable, higher order spin-spin interactions^ 1 ^, 
whereas the magnetic field applied in the plane only modifies the on-site energies of QDs: 
for B y field, QDs 1 and 3 energy level as E m {B) = E m {0) + l/2u c {R/2) 2 + gfiBS y , where 
R is the spatial separation of QDs 1 and 3 and u c = eB y /m* e c is the cyclotron frequency, 
while the QD 2 energy level only acquires the Zeeman term. Hence magnetic field in the y 
direction effectively lowers the energy of the QD 2 by AE 2 = -l/2u c (R/2) 2 with respect to 
QDs 1 and 3. 

Fig.(T5b) shows the energy spectrum of the coded qubit as a function of uj c obtained in 
LCHO-CI method. In our LCHO-CI calculation, we use the following parameters for a 
specific resonant TQD. The inter-dot distance is 10 a* B , where a* B = 9.79nm is the effec- 
tive Bohr radius for GaAs. The confining Gaussian potential on the i — th QD is of the 
form, Vi(x,y) = -V^exp (—((x - Xi) 2 + (y - yi) 2 )/d 2 ), where V® = 4.0Ry* and di = 2.5a* B . 
Ry* = 5.93777.61^ is the effective Rydberg for GaAs. In Fig.([2b), the optimal cyclotron 
frequency u* corresponds to B* = 0.69T. 

At B = 0, the four-fold degenerate ground state is separated from the four- fold degenerate 




(4) 
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spin polarized excited state by J e ff- The two computational Hilbert spaces corresponding 
to S y = ±1/2 separate energetically with increasing magnetic field while the spin polarized 
states decrease in energy. At magnetic field B*, such that J e ff = 2gfiB*, the energy gap 
is maximized and should correspond to the working point that can maintain the longest 
coherence time of the coded qubit. The coded qubit should operate at this value of magnetic 
field. However, as discussed above, the magnetic field effectively biases the QD 2. This 
removes the degeneracy of the two qubit levels and rotates them from their zero magnetic 
field states. The energy splitting of the two coded qubit levels as a function of applied 
magnetic field is shown in Fig. The splitting is a fraction of the large energy scale J e ff- 
For the largest gap, i.e. B = B*, the splitting is A. In order to restore the degeneracy of 
the two coded qubit levels, one can apply voltage to the QD 2. In Fig. 03b), microscopic 
calculations done in LCHO-CI show that a positive voltage bias on the QD 2 can indeed 
bring the TQD back on resonance in the presence of B field. We have now established the 
coded qubit and the best conditions for its operation. 

In order to operate the coded qubit, we need to be able to initialize it. We propose 
to initialize the coded qubit by turning off both interaction between QDs 1 and 3 and 
interaction between QDs 1 and 2. The only remaining interaction is between QDs 2 and 3. 
The ground state of a TQD, \L ) = | li)]^), becomes a product of a spin down state of an 
electron on the QD 1 and a singlet state of electrons across QDs 2 and 3. This intuitively is 
a ground state in magnetic field. The singlet state of two electrons in a pair of QDs 2 and 3 
can be generated in real time starting from two electrons in a biased QD 3. This procedure 
does not generate directly the coded qubit levels \q + ) and \q~), but the state \Lq) is a linear 
superposition of the two qubit levels: \L ) = ^|(e 427r / 3 |g_) — e^ l2w ^ 3 \q + )). Once interactions 
are turned on, any state can be obtained from the initial state. This can be seen by writing 
the TQD Heisenberg Hamiltonian in the basis of the two coded qubit levels {|c?±)} , 



If we take J i3 = J 2 3, and let 2 J 12 > J13 + J23, then the Heisenberg Hamiltonian corresponds 
to o x operation. If we take J23 > J13 and 2J 12 = J13 + J23, then the coded qubit Hamiltonian 
Eq.© corresponds to b y operation. The capability to rotate a qubit with respect to two 
different axes on a Bloch sphere allowes us to generate arbitrary single qubit operations. 
In practice, we tune the exchange interaction through biasing QDs and changing their 
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energies. As long as the biasing \AE\ « U — V is satisfied, we expect the quantum state 
to remain in the qubit subspace during the process of tuning the exchange interactions as 



discussed in Ref. 
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Next, we discuss the measurement of coded qubits. Several proposals^^yi discusses 
methods of detecting chirality of a triangular (three-body) antiferromagnetic cluster. For 
our specific TQD-based coded qubit, we propose to apply a B z field to split the two coded 
qubit levels with different chirality. As discussed earlier, the application of B z field will 
give rise to additional terms proportional to chirality operator on top of the Heisenberg 
model, Eq.([2J), used for qubit modeling. Since the additional term to the Heisenberg model 
is strictly proportional to the chirality, it splits the two coded qubit levels with the gap given 
by 



(u - vy- v *c 

where $ is the magnetic flux through the device and $o = eh/c is the flux quanta. This en- 
ergy gap can be interpreted as the energy difference between two magnetic dipole moments^ 
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Scarola et. al. also showed that 



oriented in opposite directions under B z . In Ref. 
the unintended qubit rotation due to the additional chirality term can be significantly re- 
duced for a range of specific weak B z fields that satisfy R/lu > 1 and u c /u < 1, where 
l h = ^Jhc/eBil + 4u^/o;2)-i/4 is the mo dified magnetic length, and u a is the confining 
frequency used to approximate the potential of a QD. Thus, initializing the system under 
B z might not be a good idea as the field has to stay on for the entire period of quantum 
computation and eventually modify the quantum state of the coded qubit due to magnetic 
moment coupling with B z . However, a measurement done with advanced spectroscopy only 
requires B z to be turned on for a relatively short time and significantly limit the extent to 
which the coded qubit will be modified. 



IV. DOUBLE QUBIT OPERATIONS 

We now turn to discuss double coded qubit operations using both Hubbard and Heisen- 
berg models. Starting with the Hubbard Hamiltonian, Eq.(TI|), we derive perturbatively a 
Heisenberg Hamiltonian for a complex of six electron spins with a set of exchange interactions 



S 



Jij 



Jl3 — J23 



J12 


— J56 — 


J a = 


J45 


— J46 — 


Jb = 




J34 — 


Jc = 



At 2 

U -V 

2t 2 2t 2 



(U-V + V) u-v-v 
At' 2 



(7) 



U-V 

with all other exchange interactions set to zero. For all the Hubbard parameters considered 
in this study (i.e. weak inter-TQD tunneling and strong Coulomb interaction), J c is usually 
about two order of magnitudes smaller than J a and J&, mainly due to the fact that if < t/2 
in our study. 

We treat the inter-TQD exchange interaction, scaled by J c , perturba- 
tively to derive an effective interaction in the coupled coded qubits subspace, 
{i<Z+)i<7+}> \q+)\q-), \q!l)\q+), where the superscript L and R stand for left 

and right TQD respectively. The effective qubit-qubit interacting Hamiltonian reads 

H 2q = A\a'\(n ■ 5i)(n • S 2 ) + 2\(5'\m ■ S x + 2\f3'\m ■ S 2 , (8) 



where a' = e~^~ J c /9, and 0' = {—J a /2 + J b /2 + J c /27) , and the direction vectors in 
Eq.® are defined as follows: n = (^y/\a'\ + Re(a') / y/2\a'\, y/\oc'\ — Re(a')\^2\a'\, oV rh = 
(Re(P')/\j3'\, Im(f3')/\(3'\, 0). Thus, the effective interaction for 2 TQD-based code qubits is 
equivalent to an XY Hamiltonian under a uniform in-plane magnetic field for 2 spins. 

If we rotate the effective qubit-qubit interacting Hamiltonian, Eq.(jS]), from the coded 
qubit basis to the Jacobian basis, {\Lq)\Lq), |Lq)|L^ ), \Li)\Lq), \Li)\Lf )}, the rotated 
Hamiltonian reads 

H 2Lq = aSlS z 2 + 0S1 + (9) 

where a = J c /9, and (3 = —J a /2 + Jb/2 + J c /18 . The Jacobian state \L ) is defined in 
Sec JIIII when we discuss the initialization of the coded qubit. The Jacobian state \L±) = 



a/1/3| 4,1) ^23) — a/2/3| ti) l^23)> where |T 2 ° 3 ) is a S z — triplet state of electrons on QDs 
2 and 3 and \T 2 ~ 3 ) is a S z = —1 triplet state of electrons on QDs 2 and 3. Similarly, \Li) = 
-l/V2(e i2w / 3 \q^) + e- l2n / 3 \q+)) can also be written as a linear combination of the coded 
qubit levels. Therefore this rotation of basis only requires single qubit operations applied 
to each coded qubit; the two Hamiltonians, Eq.® and Eq.©, are locally equivalent 35 . As 
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Ising interaction can be used to generate CNOT gate^, we show how to generate non-local 
two qubit interactions with 2 TQD-based coded qubits. 

Fig. 01]) shows results of exact diagonalization of the Hubbard Hamiltonian for two coupled 
TQDs as a function of increasing coupling t' between the two molecules. The inset shows 
the entire energy spectrum at t'/t = 0.2. The low energy spectrum consists of four levels, 
characterized mostly by the two lowest levels of each TQD. These 4 low-lying energy levels 
are well separated from the rest of the spectrum and constitute the two coupled coded qubit 
subspace. The low-lying spectrum contains a doubly degenerate level, which is a signature 
of the Ising model in an external field. In practice, the tuning of J c , the inter-TQD exchange 
interaction, should be done via tuning the tunneling parameter t'. Fig. fll]) shows that the 
Ising model features of the energy spectrum are maintained over a wide range of values of t', 
and hence the Ising model Hamiltonian, derived above, describes the coupling of two coded 
qubits very well. 



V. CONCLUSION 



In summary, we present a theory of quantum circuits based on coded qubits encoded 
in chirality of electron spin complexes in lateral gated semiconductor triple quantum dot 
molecules with one electron spin in each dot. Using microscopic Hamiltonian and exact 
diagonalization techniques we show how to initialize, coherently control and measure the 
quantum state of a chirality based coded qubit using static in-plane magnetic field and 
voltage tuning of individual QDs. The microscopic model of two interacting coded qubits is 
established and mapped to an Ising Hamiltonian. Hence both conditional two-qubit phase 
gate and voltage controlled single qubit operations are demonstrated. 
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(a) 




(b) 





(c) 




FIG. 1: Three possible layouts for a chain of TQD-based coded qubits. Each solid and dash line 
represents an exchange interaction between a pair of QDs. The layouts (a) and (b) differ in that 
the constituent coded qubits are implemented with triangular TQD in (a) and with linear TQD 
in (b). The layouts (a) and (c) differ in the way logical qubits are connected in a chain. If the 
inter-qubit distance is reduced beyond a threshold, an additional set of exchange interactions will 
arise and dominate the behavior of the system in layout (c). 
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FIG. 2: (a). A schematic of a single TQD initialized in S y = — 1/2 subspace under a B y field. The 
By field breaks the discrete rotational symmetry, and it affects QD 2 differently from QDs 1 and 3. 
(b). The energy spectrum of a triangular TQD as a function of lj c /uj*, where u c is the cyclotron 
frequency used in LCHO-CI calculation, uj* corresponds to a magnetic field B* such that the gap 
between computational subspace (the lowest e|i4rgy level in the plot) and the rest of spectrum is 
maximized as indicated by the black arrow. 
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FIG. 3: (a). The energy spectrum of qubit levels (the zoom in of the lowest level in Fig. (2b)) as 
a function of u c /oj*. The B y field splits the doubly degenerate ground state, and A corresponds to 
the gap between the two levels at co* or equivalently at B*. (b). The energy spectrum of the qubit 
levels under B* as a function of voltage V2 on QD2, and V^ i s the unbiased voltage at resonant 
condition when there is no external B field. The splitting of levels under B field can be restored 
via gate voltage tuning. 
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FIG. 4: The four lowest energy levels of 2 coupled triangular TQDs as a function of inter-TQD 
tunneling t' . The energy spectrum resembles that of the Ising model with an external field: a 
doubly degenerate levels corresponding to state | |t) an d | tl) an d two unique levels | H) and 
| ft). "X2" denotes a doubly degenerate level. Inset: The entire energy spectrum at t'/t = 0.2 
calculated with LCHO-CI method. 
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